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We continue our lattice simulations of QCD with 2 flavours of colour-sextet quarks 
as a model for conformal or walking technicolor. A 2-loop perturbative calculation 
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of the /3-function which describes the evolution of this theory's running coupling 

CO ! constant predicts that it has a second zero at a finite coupling. This non-trivial zero 

> 

ON ' would be an infrared stable fixed point, in which case the theory with massless quarks 

m 



would be a conformal field theory. However, if the interaction between quarks and 



. antiquarks becomes strong enough that a chiral condensate forms before this IR fixed 

o 

point is reached, the theory is QCD-like with spontaneously broken chiral symmetry 
and confinement. However, the presence of the nearby IR fixed point means that 
there is a range of couplings for which the running coupling evolves very slowly, i.e. it 
'walks'. We are simulating the lattice version of this theory with staggered quarks at 
finite temperature studying the changes in couplings at the deconfinement and chiral- 
symmetry restoring transitions as the temporal extent (Nt) of the lattice, measured 
in lattice units, is increased. Our earlier results on lattices with iVj = 4, 6 show 
both transitions move to weaker couplings as Nt increases consistent with walking 
behaviour. In this paper we extend these calculations to Nt = 8. Although both 
transition again move to weaker couplings the change in the coupling at the chiral 
transition from Nt = 6 to Nt = 8 is appreciably smaller than that from Nt = 4 to 
Nt = 6. This indicates that at Nt = 4, 6 we are seeing strong coupling effects and 
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that we will need results from Nt > 8 to determine if the chiral-transition coupling 
approaches zero as Nt — > oo, as needed for the theory to walk. 

I. INTRODUCTION 

We are interested in extensions of the standard model which have a strongly-coupled 
(composite) Higgs sector. The most promising theories of this type are the so-called tech- 
nicolor theories Q] , QCD-like gauge theories with massless (techni-)quarks, where the 
(techni-)pions play the role of the Higgs field, giving masses to the W and Z. Such theories 
tend to have phenomenological problems, especially when they are extended to give masses 
to quarks and leptons. Walking technicolor theories, where gauge group and fermion con- 



tent are chosen so that the runnin 
be able to avoid such difficulties [3 



coupling constant evolves very slowly ('walks'), might 
6|. Deciding whether a candidate gauge theory has the 
properties needed is a non-perturbative question. Hence lattice gauge theory simulation 
methods are the only way to answer this reliably. 

For a given gauge group with Nf fermions in a specified representation of that group, 
there is some value of Nf, below which the gauge theory is asymptotically free. Below this 
value there is a range of Nf for which the second term in the perturbative Callan-Symanzik 
/3-function has the opposite sign from the first. Hence, if the 2-loop f3 function describes 
the physics the theory with Nf in this range, ft has a second non-trivial zero representing 
an infrared (IR) fixed point. If this is true the theory is a conformal field theory with a 
continuous spectrum. However, there is a second possibility. If the fermion-antifermion 
coupling becomes strong enough that a chiral condensate forms before the would-be fixed 
point is reached, this effectively removes the fermions from consideration for longer distances, 
the IR zero is avoided, and the coupling approaches infinity at large distances. In this case 
the theory is QCD-like with confinement as well as chiral-symmetry breaking. However, the 
presence of the nearby IR fixed point means that the /3 function becomes small at some 
value of the coupling, and the coupling constant 'walks'. 

If we restrict ourselves to SU(N C ) gauge groups with N c relatively small, there are a 
limited number of potential candidates. These have been identified and rough estimates of 

D 

the value of Nf, which separates conformal from walking behaviour, have been made [7h 
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Extensive lattice studies have been made for N c = 3 with fermions in the fundamental 
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representation of the colour group. [15l436| There have also been studies with N c = 2 and 
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as well as studies with 



fermions in the fundamental representation of colour 
N c — 2 and fermions in the adjoint (symmetric tensor) representation 39N48j| . Finally, there 
have been studies with N r = 3 and fermions in the sextet (symmetric-tensor) representation 
of the gauge group. [49M56] . 

We are concentrating our efforts on QCD (iV c = 3) with colour-sextet quarks. For this 
choice, asymptotic freedom is lost at Nf = 3^. This means that only Nf = 2,3 are of 
interest. Both of these have /3- functions where the 1 and 2-loop terms are of opposite sign. 
Nf = 3 is close enough to the number of flavours for which asymptotic freedom is lost that 
the IR fixed point occurs for very weak coupling, for which perturbation theory can probably 
be trusted. Hence it is believed that this theory is almost certainly a conformal field theory. 
Estimates of the value of Nf, which separates conformal from walking behaviour, suggest 
that Nf = 2 is a good candidate for walking behaviour. However, there is enough uncertainty 
in these methods that a more reliable study of the Nf = 2 is waranted. We have thus chosen 
to study the Nf = 2 theory, using lattice gauge simulations with staggered quarks. Lattice 
studies of QCD with 2 flavours of colour-sextet Wilson quarks have been performed by 
Degrand, Shamir and Svetitsky. To date these have been unable to tell unambiguously 
whether this theory is conformal or walking. The Lattice Higgs Collaboration have been 
studying this theory using improved staggered quarks. Recently they have reported evidence 
that this theory spontaneously breaks chiral symmetry which would indicate that it has 
walking behaviour (unless there is a bulk chiral transition at even weaker coupling). 

Whereas the other groups have concentrated their efforts on determining the nature of 
QCD with 2 sextet quarks from studies of the zero temperature behaviour of the theory 
(apart from some early thermodynamics simulations by Degrand, Shamir and Svetitsky), 
we are studying the thermodynamics of this theory. Here we are measuring the depen- 
dence of the lattice (bare) coupling at the deconfinement and chiral-symmetry restoration 
transitions on N t , the temporal extent of the lattice in lattice units. If these are indeed 
finite-temperature transitions, the couplings at which they occur should tend to zero as 
N t — y oo in a manner controlled by asymptotic freedom. Such behaviour would indicate 
walking. If the theory is conformal, these couplings should approach a non-zero constant as 
Nt — > oo, indicating a bulk transition. Simulations at N t = 4 and 6, reported in our earlier 
publication showed that both transition couplings did decrease with increasing N t . This 
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work reports the results of simulations at N t = 8. While both transitions do tend to weaker 
couplings as N t goes from 6 to 8, the change in coupling at the chiral transition, which occurs 
at a considerably weaker coupling that the deconfmement transition, is much smaller than 
that between N t = 4 and 6. (Such separation of the deconfmement and chiral-symmetry 
restoration transitions, which is not observed for fundamental quarks, has been observed 
with adjoint quarks 57|, |58|). The most likely interpretation is that between N t = 4 and 6, 
this transition is in the strong-coupling domain where the quarks have condensed to form a 
chiral condensate at length scales of order of the lattice spacing and do not participate in 
the running of the coupling constant, which now runs as in quenched QCD. Between N t = 6 
and 8 the chiral transition coupling finally emerges into the weak-coupling regime where the 
quarks also particpate in the running of the coupling constant. This means that we will need 
to simulate at even larger N t s to determine whether this theory is conformal or walking. 

As for Nt = 6, the N t = 8 lattice shows a clear 3-state signal above the deconfmement 
transition. These states are characterized by the phase of the Wilson Line (Polyakov Loop) 
having the values 0, ±^f-, a vestige of the Z3 colour symmetry of the quenched theory. Within 
the limitations of our simulations, all 3 states appear stable. At even weaker couplings - 
close to the chiral transition - the 2 states with complex phases disorder to a phase with 
a negative Wilson Line. This phase structure, which is richer than that for fundamental 
quarks, where the Wilson Line is always real and positive, was predicted by Machtey and 
Svetitsky and observed in their simulations with Wilson quarks 59]. 

In section 2 we describe our simulation techniques, and how one can measure the running 
of the coupling constant from thermodynamics. Section 3 describes our simulations and 
results. Finally in section 4 we discuss our results, draw conclusions, and indicate directions 
for future investigations. 



II. METHODOLOGY 



For the gauge fields we use the standard Wilson (plaquette) action: 
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For the fermions we use the unimproved staggered- quark action: 
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where Tp = £^ Vii D n with 

D^(x) = l -[Uf\x)^{x + p)- U^(x - fiMz - p.)], (3) 

where is the sextet representation of U, i.e. the symmetric part of the tensor product 
U ®U . When Nf is not a multiple of 4 we use the fermion action: 

Sf = E X j {[P+m][-p+m]} N f/ 8 X . (4) 

sites 

The operator which is raised to a fractional power is positive definite and we choose the real 
positive root. This yields a well-defined operator. We assume that this defines a sensible field 
;heory in the zero lattice-spacing limit, ignoring the rooting controversy. (See for example 



60[ for a review and guide to the literature on rooting. 



We use the RHMC method for our simulations [6l|, where the required powers of the 
quadratic Dirac operator are replaced by diagonal rational approximations, to the desired 
precision. By applying a global Metropolis accept/reject step at the end of each trajectory, 
errors due to the discretization of molecular-dynamics time are removed. 

Finite temperature simulations are performed by using a lattice of finite extent N t in 
lattice units in the Euclidean time direction, and of infinite extent N s in the spatial direction. 
In practice this means we choose N s ^> N t . The temperature T = 1/N t a, where a is the 
lattice spacing. (In our earlier equations we set a = 1.) Since the deconfinement temperature 
Td and the chiral symmetry restoration temperature T x should not depend on a, and since 
a = 1/NfT, measuring the coupling g at T d or T x as a function of N t gives g(a) for a series 
of a values which approach zero as N t — > oo. If the ultraviolet behaviour of the theory is 
governed by asymptotic freedom, g(a) should approach zero as a — > 0, i.e. N t — » oo. The 
way gd and g x approach zero should be determined by the perturbative (3 function. The 
2-loop (3 function 

0(g) = -b l9 3 - b 2 g 5 . (5) 

Then expressing our coupling constant evolution in terms of /3 = Q/g 2 (We apologize for the 
fact that we are using (3 for 2 different purposes) 

AP(P) = 0(a) - 0(\a) = (12&i + 72b 2 /0) ln(A) (6) 

through this order. For Nf flavours of sextet quarks, 

10„\ ... 2 
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102 - 2 -f Nf j /(16vr 2 ) 2 (7) 

If, on the other hand, the Nf = 2 theory is conformal, the continuum, zero coupling 
{(3 — > oo) limit has an unbroken chiral symmetry (and is unconfined). Hence there will be 
a bulk chiral transition at a finite coupling, which survives in the N t — > oo limit, so the 
coupling and hence (3 at the chiral transition will tend to a finite value in this limit. (Since 
the (3 value at the deconfinement transition (/3d) is expected to be less than that at the chiral 
transition (f3 x ), it follows that [3d will also approach a finite value as Nf — > oo.) 

We determine the position of the deconfinement transition as that value of (3 where the 
magnitude of the triplet Wilson Line (Polyakov Loop) increases rapidly from a very small 
value as [3 increases. The chiral phase transition is at that value of (3 beyond which the chiral 
condensate (ipip) vanishes in the chiral limit. Because we are forced to simulate at finite 
quark mass, this value is difficult to determine directly. We therefore estimate the position of 
the chiral transition by determining the position of the peak in the chiral susceptibility Xm 
as a function of quark mass, and extrapolating to zero quark mass. The chiral susceptibility 
is given by 

= V [<(^) 2 > - (W 2 ] (8) 

where the () indicates an average over the ensemble of gauge configurations and V is the 
space-time volume of the lattice. Since the fermion functional integrals have already been 
performed at this stage, this quantity is actually the disconnected part of the chiral sus- 
ceptibility. Since we use stochastic estimators for ipip, we obtain an unbiased estimator for 
this quantity by using several independent estimates for each configuration (5, in fact). Our 
estimate of (ipip) 2 is then given by the average of the (10) estimates which are 'off diagonal' 
in the noise. 

Our N t = 8 simulations are performed on 16 3 x 8 lattices. Near the chiral transition, 
where finite size effects are a concern, we also perform simulations on a 24 3 x 8 lattice for 
the lowest quark mass. We perform simulations with quark masses m = 0.005, m = 0.01 
and m = 0.02 in lattice units, to enable continuation to the chiral (m = 0) limit. (Since 
we do not have any zero temperature measurements, the more desirable method of choosing 
lines of constant physics is impossible.) Our trajectory length is chosen to be At = 1 where 



t is the molecular-dynamics 'time' in HEMCGC normalization 



62). 



A more detailed discussion of our methods of choosing parameters, run lengths, etc. is 



given in our earlier paper describing our N t = 4, 6 simulations 



53]. 



III. SIMULATIONS AND RESULTS 

We simulate QCD with 2-flavours of colour-sextet staggered quarks on 16 3 x 8 and 24 3 x 8 
lattices. For the smaller lattice we perform simulations with masses m = 0.005, m = 0.01 
and m = 0.02 to allow extrapolation to the chiral limit, for a set of (3 values covering the 
range 5.5 < (3 < 7.4. To probe the various phases of the Wilson Line, we use 2 different sets 
of runs. In the first set of runs we use an ordered start, in which the gauge fields are set to 
the unit matrix on all links, at the highest (3, and use configurations from higher /3s to start 
runs at lower /3s. The second set of runs uses a start in which the gauge fields are set to the 
unit matrix, except for the timelike gauge fields on a single timeslice, which are set to the 
matrix diag(l, —1, —1). This puts the system in a state with a real negative Wilson Loop 
at large /3 s. 

The length of a typical run at a fixed (/3,m) away from the transitions is 10,000 trajec- 
tories. Close to the deconfinement transition, this is increased to 50,000 trajectories. Run 
lengths of 50,000 trajectories are also used close to the transition from a state where the 
Wilson Line has phase ±27r/3 to one where it has phase tt. We have detailed our run lengths 
in the appendix. 

Since finite (spatial) volume effects are most likely to be present in the weak-coupling 
domain at small quark masses, where they have the potential to shift the chiral transition, 
we have also performed a set of simulations on 24 3 x 8 lattices at the lowest quark mass. 
These simulations at m = 0.005 cover the range 6.2 < /3 < 7.4 with mesh 5(3 = 0.1, and 
with 10,000 trajectories at each /3, from positive Wilson line starts. 



A. Results 



Starting from large (3 values, the runs which start from a completely ordered state with 
Wilson Line +3 continue to have positive Wilson Lines down to (3 = 5.8 for m = 0.02, 
and down to (3 = 5.7 for m = 0.01 and m = 0.005. Below these (3 values, which are just 
above the deconfinement transition, we see a clear 3-state signal, where the system tunnels 
between states where the Wilson Line has phases 0, ±27r/3. Because of this, we bin our 
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data according to the phase of the Wilson Line for each configuration. Configurations 
where — 7r/3 < < 7r/3 are considered to be in the = bin. Outside of this range the 
configurations are considered to be in the ±27r/3 bins depending on whether the imaginary 
part of the Polyakov loop is positive or negative. These last 2 bins are combined by complex 
conjugating those Wilson Lines which have negative imaginary parts. 

Starting from large /3s, in those runs which start from the second ordering with Wilson 
Line —1, the Wilson Line remains negative down to /3 ~ 6.9 for m = 0.02, « 6.8, m = 0.01 
and (3 rs 6.7, m = 0.005. Below these values the system makes a transition to a state with 
Wilson Line phase ±27r/3. Below these /3 values, these runs remain in states with Wilson 
Line phases ±27r/3 down to (5 = 5.8, for each m. For /3 = 5.7 and below we see clear 3-state 
signals where the system tunnels between the 3 states. For this reason we again bin our 
data according to the Wilson Line phase, 0. 

In figure we present the Wilson Line and chiral condensate for the states with a real 
positive Wilson Line plotted against (3, for each of the 3 masses. In figure [lb we plot the 
magnitude of the Wilson Line and the chiral condensate for those states with complex or 
real negative Wilson Lines. 

In both sets of graphs in figure [1] we observe a rapid increase in the (magnitude of the) 
Wilson Line at ~ 5.65, corresponding to the deconfinement transition. The sudden drop 
in the magnitudes of the Wilson Lines of figure [lb for 6.6 < /3 ^ 7.0 marks the transition 
where states with complex Wilson Lines (0 = ±27r/3) disorder to a state with a real negative 
Wilson Line = it. 

It is clear that the chiral condensate becomes small for large 0, and decreases with 
decreasing quark mass, which suggests that it will vanish in the chiral limit, for /3 large 
enough. However, extrapolating the chiral condensate to zero quark mass to determine 
the chiral transition from these quark masses where the (3 dependence is so smooth would 
be exceedingly difficult. We therefore estimate the position of the chiral transition from 
determinations of the positions of the peaks in the chiral susceptibilities for each mass. These 
are plotted in figure Eh- For the lower two masses, the peaks in the chiral susceptibilities 
are well defined. (This is the best evidence we have that our quark masses are small enough 
to perform the chiral extrapolation.) In addition, for the limited set of /3 values of our 
simulations, both the m = 0.01 and the m = 0.005 'data' peak at the same /3, namely 
(3 = 6.7. We therefore estimate that the position of the peak and thus the chiral phase 
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FIG. 1: a) Wilson Line and chiral condensate for real positive Wilson Line states as functions of j3 
on a 16 3 x 8 lattice. 

b) Magnitude of Wilson Line and chiral condensate for states with complex or real negative Wilson 
Lines as functions of /3 on a 16 3 x 8 lattice. 
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transition at m = are at X = 6.7(1). This means that (3d and (3 X are well separated as 
was observed for N t = 4 and 6. At m — 0.005, close to the chiral transition, we have also 
performed simulations on larger (24 3 x 8) lattices. The Wilson Lines and chiral condensates 
show little difference between the two lattice sizes. The chiral susceptibilities plotted in 
figure [2] indicate that finite size effects are indeed small. This is more significant, since such 
fluctuation quantities are most sensitive to finite volume effects. 

We have also looked at the chiral susceptibilities for the states with real negative or 
complex Wilson lines and find peaks at (3 = 6.8(1) for m = 0.02, f3 = 6.7(1) for m = 0.01 
and (3 = 6.6(1) for m = 0.005. Since these values are close to the transitions from the state 
with a negative Wilson Line to states with complex Wilson Lines, there is a possibility of 
interference between these two transitions. For this reason we concentrate our studies on 
the chiral transition measured in the positive Wilson Line state. 

We now turn our attention to more precise estimates of (3 d . For this purpose, we histogram 
the magnitudes of the Wilson Lines in the neighbourhood of the deconfinement transition. 
Such histograms are shown in figure [3] for each quark mass. We estimate that the transition 
occurs aX (3 = (3d = 5.66(1) for m = 0.02, at {3d = 5.65(1) for m = 0.01 and at fid = 5.65(1) 
for m = 0.005. 

The positions of the deconfinement and chiral transitions, ex trap olated to the chiral limit 
are given in tabled for each of the 3 N t values (N t = 4, 6 from 53j, N t = 8 this work). 



N t 


Pd 


Px 


4 


5.40(1) 


6.3(1) 


6 


5.54(1) 


6.6(1) 


8 


5.65(1) 


6.7(1) 



TABLE I: Nf = 2 deconfinement and chiral transitions for Nt = 4, 6, 8. 



B. Interpretation 



We now compare the changes in /3 d and (3 X with what would be expected if the running of 
the lattice coupling constant is given by the 2-loop (perturbative) (3 function - equation [6j 
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FIG. 2: a) The chiral susceptibilities as functions of j3 for each of the 3 masses on a 16 3 x 8 lattice, 
b) The chiral susceptibilities at m = 0.005 as functions of j3 on a 16 3 x 8 lattice and on a 24 3 x 8 
lattice. 



12 



16°x8 lattice m=0.005 



16 J x8 lattice m=0.01 



800 



,9=5.64 0=5.65 ,3 = 5.66 /S=5.67 




0.00 0.05 0.10 

|WILS0N LINE | 

800 



0.15 



0.00 0.05 0.10 

|WILSON LINE | 



16°x8 lattice m=0.02 




0.05 0.10 
|WILSON LINE| 



0.15 



FIG. 3: Histograms of magnitudes of Wilson Lines: a) for m = 0.005, b) for m = 0.01 and c) for 
m = 0.02. 

For the deconfinement transition 

A&(6,4) = d (N t = 6) - p d (N t = 4) « 0.14 (9) 

compared with the prediction of equation [6] which predicts A$j(6, 4) w 0.12, whereas 

A/3 d (8, 6) = /3 d (JV t = 8) - ^(JV f = 6) « 0.11 . (10) 



compared with the 2-loop prediction A/3 d (8, 6) ~ 0.09. If the quarks were actively screening 
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colour at these couplings (5.40 < (3 < 5.65), it would not be unreasonable to assume 
that these deconfmement couplings were weak enough to be governed by the perturbative 
P function. The fact that the measured A/3s are within ps 20% of those predicted by 2- 
loop perturbation theory would tend to support this interpretation. However, examining 
the running of the coupling constant at the chiral transition, will lead us to a different 
conclusion. 

For the chiral transition, we find 

A/3 X (6,4) = P x (N t = 6) - P x (N t = 4) ps 0.3 (11) 

while 

A/3 X (8, 6) = /3 x (N t = 8) - P x (N t = 6) ps 0.1 . (12) 
Using equation [6] to estimate A/3 X , yields 

A/3 X (6,4) ps 0.122 (13) 

and 

A/3 X (8,6) ps 0.087. (14) 

While A/3 X (8,6) appears consistent with our measurements, A/3 X (6, 4) does not. What this 
suggests is that for N t in the range 6-8, /3 X is in the weakly-coupled domain where scaling 
is controlled by asymptotic freedom, while N t in the range 4-6 is in the strongly coupled 
domain. 

In the strongly coupled domain, the fermions have formed a chiral condensate, which 
effectively stops them from contributing significantly to the running of the coupling constant. 
Hence we expect that the running of the coupling in this region will be that of the quenched 
theory, i.e. that for equation |6] with Nf = 0. This yields 

A/3 X (6,4) ps 0.357, (15) 

which is consistent with what we observe. (It also gives 

A/3 X (8,6) ps 0.253, (16) 

which is larger than what we observe.) Thus we conclude that the chiral transition emerges 
from the strongly coupled domain, where the quarks play little part in the coupling constant 
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evolution, into the weak coupling regime, where the running of the coupling is determined 
by asymptotic freedom, around (3 x (N t = 6). 

One might argue that both A/3 x s are consistent with either Nf = or Nf = 2 scaling, 
because of the relatively large error-bars in table [B However, comparing the graphs of the 
chiral condensates for fixed masses for each N t - see figure H]- we see that A/3 X (6,4) really 
does appear to be much larger than A/3 X (8,6), and that the estimates of equations 11 and 
12 are more accurate than the errors in the individual /3 x s would suggest. 
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FIG. 4: Chiral condensates for m = 0.005 for Nf = 4, 6 and 8, in the high /3 (weak coupling) 
regime. 



This interpretation of the running of the lattice coupling constant at the chiral transition 
indicates that the region /3 ^ 6.6 is one of strong coupling, governed by quenched dynamics. 
In particular, the change in (3 d as N t is varied from 4 to 6 to 8 should be governed by quenched 
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dynamics. However, it has been determined that the evolution of the deconfmement coup ling 
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64). 



in quenched QCD, is only well described by 2-loop perturbation theory for > 6.1 
Hence the changes in [3d that we observe are not expected to be described by quenched 
perturbation theory. We note, however, that (3d(N t = 6) = 5.54(1) is close to (3d(N t = 
3) ~ 5.55 found for the quenched theory 65^ , w hile (3d{N t = 8) = 5.65(1), compared to 
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671 ]. Since the ratios of lattice spacings 



0d(N t = 4) 5.69 for the quenched theory 
in these 2 cases are identical, such comparison is justified. Taking into account the fact 
that small N t effects are expected to make the N t = 3 quenched (3d anomalously small, 
the comparison is remarkably good. Unfortunately we cannot expect similar comparisons 
between f3d{N t = 4) and f3d{N t = 2) for the quenched theory, to work, since it is well 
known that Pd{N = 2) « 5.1 for quenched lattice QCD, is anomalously low [65]. Still it 
is reasonable that a strong-coupling quenched interpretation of the running of fid is correct 
for Nt = 4, 6, 8. The fact that the changes in (3d between N t = 4, 6 and 8 are appreciably 
less than predicted by quenched perturbation theory is a well-known feature of the strong- 



coupling domain of quenched lattice QCD 



63. 



64]- 



IV. DISCUSSION AND CONCLUSIONS. 



We are studying thermodynamics of QCD with 2 massless colour-sextet 'quarks' in an 
attempt to distinguish whether it is a conformal field theory or if it 'walks', i.e. is a confined, 
chiral-symmetry broken theory, with a slowly evolving coupling. Simulations are performed 
on lattices with spatial extent N s a and temporal extent N t a (a is the lattice spacing) with 
N s » N t . We use staggered quarks with several small masses to extrapolate to the massless 
quark limit. The temperature T = 1/N t a. Hence for fixed T, chosen to be either the 
deconfmement temperature Td or the chiral-symmetry restoration temperature T x , we can 
vary a by varying Nt, thus studying the running of the coupling g(a) as a — > by simulating 
at a series of N t s increasing towards infinity. Our earlier simulations were performed using 
N t = 4, 6. Those we report here use N t = 8. 

At N t = 8, as at the smaller N t s, the chiral transition occurs at a much weaker coupling 
than the deconfmement transition (see tabled]). (This contrasts with the case with quarks 
in the fundamental representation of colour, where these transitions appear coincident.) 
Between Nt = 4 and 6 both transitions move to appreciably smaller couplings. While this 
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trend continues between N t — 6 and 8, the change in couplings at the chiral transition, is 
much smaller than that between N t — 4 and 6. A possible explanation is that for couplings 
between those at the N t = 4 and the N t — 6 chiral transitions, the system is in the strong- 
coupling regime, where the quarks are bound in a chiral condensate at distances ^ a, 
and thus do not contribute significantly to the evolution of the coupling constant, which 
thus evolves as in quenched (Nf = 0) QCD. Between the couplings for the N t — 6 and 
N t — 8 transitions, the system emerges into the weak-coupling domain, where the fermions 
contribute and the coupling evolves according to Nf = 2 asymptotic freedom. Although 
we have given semi-quantitative evidence for this interpretation, we cannot rule out the 
possibility that the coupling at the chiral transition is approaching a fixed non-zero value. 
This would be evidence for a bulk chiral transition, implying that the continuum theory has 
unbroken chiral symmetry and is thus conformal. 

In order to distinguish conformal from walking behaviour, we will need to perform sim- 
ulations at larger N t values. We have already begun simulations on N t = 12 lattices. In 
addition, since the changes in /3 X between N t = 6 and N t = 8 and those expected between 
N t = 8 and Nt = 12 are of order 0.1 we will need more f3 values in the neighbourhood of 
f3 x to determine this value more accurately. We are currently performing such simulations 
and additional simulations with a smaller quark mass (m = 0.0025) at N t = 8, to aid with 
the chiral extrapolation. With these new simulations, we are concentrating on the chiral 
transition, since it would require N t values much greater than what is currently feasible to 
have (3d in the weak-coupling (f3 > 6.6) regime. 

Our runs at N t = 8 show a phase structure similar to what was observed at N t = 6. 
Above the deconfinement transition, the Wilson Line exhibits a definite 3-state structure. In 
addition to the state with a positive Wilson Line, which is all that is observed for fundamental 
quarks, there are states characterized by Wilson Lines oriented (at least approximately) in 
the directions of the 2 complex cube roots of unity. The existence of this 3-state signal is 
probably because chiral symmetry is still broken in this regime, effectively removing the 
fermions from the dynamics, so that it behaves as a quenched theory. This suggests that 
the 3-state signal is the vestage of the spontaneously broken Z 3 symmetry of the deconfined 
pure gauge theory. The fact that this Z 3 symmetry is explicitly broken manifests itself in 
the fact that the magnitude of the Wilson Line in the complex Wilson Line states is smaller 
than that in the positive Wilson Line state. Within the limitations of our simulations, all 
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3 states appear to be stable. At some (3 value close to the chiral transition, the complex 
Wilson Line states disorder to a state with a negative Wilson Line. Above this transition 
the magnitude of the Wilson Line in the negative Wilson Line state is around 1/3 of that 
for the positive Wilson line state. This leads us to speculate that the transition indicates a 
breaking of colour SU(3) to colour SU{2) x U{1)y- Arguments for the existence of states 
with Wilson Lines having phases ±27r/3 and 7r in addition to that with phase have been 
given by Machtey and Svetitsky, who showed evidence for them in their simulations with 
Wilson quarks. 

We also plan simulations to measure the zero-temperature properties of this theory. In 
the weak-coupling regime (3 > (3 X , we will check whether the theory is a conformal field 
theory or if it is in the quark-gluon plasma phase of a Q CD-like gauge theory. If the theory 
is a conformal field theory (for massless fermions) all 'hadron' masses will vanish with the 
same anomalous dimension, and the chiral condensate will also vanish with an anomalous 
dimension, in the chiral limit. Because such anomalous dimensions are determined by the 
infrared attractive fixed-point, they should be independent of /3. 

If we do not find evidence of conformal behaviour, we will check for QCD-like behaviour in 
the chirally-broken phase, and for evidence that this phase has a continuum limit controlled 
by asymptotic freedom. This will require very large lattices, since we need to choose /3 
values large enough for asymptotic freedom to control the renormalization group scaling of 
observables, while keeping (3 < (3d(< (3 X ). Here we will need to measure the masses of the 
'hadrons' to determine if our quark masses are small enough and our lattices large enough to 
observe that the 'pion' masses vanish in the chiral limit proportional to y/m, while the other 
'hadrons' remain massive. We will also need to check for evidence that the chiral condensate 
remains finite in the chiral limit. In addition we will measure f„ and study propagators of 
vector and axial vector mesons which contribute to the S parameter, as is being done by 



the Lattice Strong Dynamics Collaboration, for fundamental quarks [35|. Simulations will 
be performed at several (3 values to determine the running of the bare coupling and of some 
appropriately-defined renormalized coupling. This is necessary to check that the theory has 
the correct ultraviolet completion. 

Zero temperature simulations with sextet quarks are already being performed using im- 
proved staggered quarks by the Lattice Higgs Collaboration, who presented preliminary 
results at Lattice 2010 
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Appendix A: Run details 

Table HT1 gives the length of our 16 3 x 8 runs in length-1 trajectories, for each (3 and mass. 
Where 2 numbers are given, the first is for a series of runs which started from an ordered 
configuration at large /3, while the second is from a start which gives negative Wilson Loops 
at large j3. 
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p 


m = 


= 0.005 


m - 


= 0.01 


m - 


= 0.02 


5.5 


10,000 


10,000 


10,000 


5.55 


10,000 


10,000 


10,000 


5.6 


50,000 


50,000 


50,000 


5.64 


50,000 


+ 


50,000 


50,000 


+ 


50,000 








5.65 


50,000 


+ 


50,000 


50,000 


+ 


50,000 


50,000 


+ 


50,000 


5.66 


50,000 


+ 


50,000 


50,000 


+ 


50,000 


50,000 


+ 


50,000 


5.67 


50,000 


+ 


50,000 


50,000 


+ 


50,000 


50,000 


+ 


50,000 


5.68 


50,000 


+ 


50,000 


50,000 


+ 


50,000 


50,000 


+ 


50,000 


5.7 


50,000 


+ 


50,000 


50,000 


+ 


50,000 


50,000 


+ 


50,000 


5.8 


10,000 




10,000 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


5.9 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


6.0 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


6.1 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


6.2 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


6.3 


10,000 




10,000 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


6.4 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


6.5 


10,000 


+ 


50,000 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


6.6 


20,000 


+ 


50,000 


10,000 


+ 


50,000 


10,000 


+ 


10,000 


6.7 


20 000 


4- 

1 


50 000 


10 000 


1 


50 000 


10 000 


4- 

1 


10 000 


6.8 


20,000 


+ 


50,000 


10,000 


+ 


50,000 


10,000 


+ 


50,000 


6.9 


10,000 


+ 


20,000 


10,000 


+ 


50,000 


10,000 


+ 


50,000 


7.0 


10,000 


+ 


10,000 


10,000 


+ 


50,000 


10,000 


+ 


50,000 


7.1 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


10,000 


+ 


30,000 


7.2 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


7.3 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


7.4 


10,000 


+ 


10,000 


10,000 


+ 


10,000 


10,000 


+ 


10,000 



TABLE II: Numbers of trajectories for each /3 and m for runs on 16 3 x 8 lattices. The first number 
is for simulations starting with positive Wilson Lines; the second is for simulations starting with 
negative or complex Wilson Lines. 
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